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CHARACTERIZATION  OF  THE  INITIAL  CLOUD  OF  CHLORINE 
RESULTING  FROM  A  90-TON  RAILCAR  INCIDENT 


1.0  INTRODUCTION 

A  release  from  a  railcar  containing  chlorine  or  other  toxic  industrial  chemical  while  transiting  near  a  city 
is  of  great  concern.  There  is  currently  disagreement  between  model  predictions  of  the  distance  that  the 
resulting  toxic  cloud  of  chlorine  vapor  travels  at  human  exposure  levels  resulting  in  death  or  requiring 
medical  treatment1  and  records  of  casualties  following  recent  incidents.2  This  paper  provides  a 
mathematical  analysis  of  the  physical  processes  occurring  following  breach  of  a  chlorine  railcar  through 
intentional  or  accidental  means.  The  relevant  toxic  cloud  formation  processes  following  creation  of  a  hole 
or  rupture  in  the  liquid  space  include  tank  equilibrium,  jet  velocity  and  flow  rate,  jet  adiabatic  expansion, 
air  entrainment  into  the  jet,  and  cloud  equilibrium.  At  this  point,  modem  dense  gas  hazard  assessment 
models  assume  the  jet  is  carried  downwind  with  the  passing  air  and  continue  to  entrain  air  and  its 
associated  mass,  momentum,  and  energy.  In  some  models,  dense  gas  slumping  calculations  precede  the 
interaction  with  the  wind.3 

This  paper  proposes  an  additional  process  that  occurs  when  the  rate  of  chemical  flow  out  of  the  tank 
exceeds  the  capacity  of  the  passing  air  to  carry  it  downwind.  The  chemical  accumulates  at  the  release 
location  and  forms  a  large,  flat  cloud  of  very  dense  aerosol  and  vapor  called  a  mist  pool.  The  mist  pool 
interacts  with  the  air  above  through  entrainment  of  chemical  into  the  air  across  the  cloud  top  area.  The 
mist  pool  behavior  is  compared  to  the  normal  jet  behavior  in  Figure  1. 


Figure  1.  Cloud  Formation  for  (a)  normal  jet  taken  up  by  the  wind,  and  (b)  mist  pool.* 

*  The  dashed  lines  show  the  progression  of  each  type  of  release  with  time. 

Normal  transport  and  diffusion  algorithms  still  apply  to  this  vapor  area  flux.  The  result  of  mist  pool 
formation  is  thus  an  area  source  of  longer  duration  than  the  tank  empty  time  leading  to  lower 
concentrations  downwind.  A  criterion  is  proposed  to  determine  when  a  mist  pool  will  form.  Attacks  on 
railcars  filled  with  chlorine  and  ammonia  are  used  as  examples  to  demonstrate  the  complete  sequence  of 
events. 
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2.0  THEORY 


2.1  Attack  and  Tank  Conditions 

The  processes  addressed  by  the  mist  pool  are  relevant  to  industrial  chemicals  that  are  stored  in  multi-ton 
quantities  as  liquid  under  pressure  and  at  temperatures  above  the  chemical  boiling  point.  Relevant 
chemicals  include  chlorine,  anhydrous  ammonia,  liquefied  natural  gas,  Freon,  propane,  and  butane.  Steam 
stored  at  temperatures  below  about  180  °C  exhibits  similar  physical  processes.  The  hole  or  rupture 
resulting  from  the  attack  is  assumed  to  be  well  below  the  chlorine  vapor-liquid  interface,  most  likely 
about  !4  of  the  way  up  from  the  tank  bottom  to  the  top.  Breaches  in  the  vapor  space  are  not  addressed. 

This  paper  will  focus  on  chlorine  railcar  tanks,  which  are  only  transported  full  or  empty.  Ammonia  stored 
in  the  same  tank  under  similar  conditions  is  considered  for  comparison  purposes. 

The  tank  cars  are  not  refrigerated,  but  the  4-inch  layer  of  insulation  between  the  tank  and  the  metal  jacket 
results  in  a  very  slow  temperature  increase  from  the  filling  temperature  to  ambient  temperature.  The 
derivations  and  analysis  here  assume  there  has  been  sufficient  time  (e.g.,  several  days)  between  filling  and 
incident  occurrence  for  the  chlorine  liquid  to  reach  ambient  temperature.  The  difference  between  the  tank 
temperature  and  the  liquid  boiling  temperature  is  defined  as  the  temperature  of  superheating.  The  range  of 
superheating  considered  is  taken  to  be  between  20  °C  and  80  °C.  Chlorine  has  a  boiling  point  of  -34  °C, 
and  ammonia  has  a  boiling  point  of  -33  °C,  so  the  relevant  tank  temperature  range  is  between  -13  °C  and 
46  °C.  The  characteristics  of  the  jet  passing  through  the  hole  change  outside  of  this  range.  Lower 
superheating  has  been  shown  to  result  in  significant  rainout  and  pooling,  while  higher  superheating  can 
result  in  vapor  forming  while  the  liquid  passes  through  the  hole  instead  of  after  exposure  to  ambient 
pressure. 

2.2  Tank  Emptying  Equilibrium 

Although  the  formation  of  jet  flow  through  the  hole  may  be  the  first  physical  process  to  occur  upon  hole 
formation,  the  jet  flow  rate  is  a  function  of  tank  pressure  and  jet  density,  both  of  which  vary  as  the  tank 
empties.  The  following  assumptions  are  made:  1)  negligible  heat  transfer  through  the  tank  and  insulation 
layers  during  the  tank  emptying  process;  2)  no  air  enters  the  tank  until  after  it  depressurizes  to  ambient 
pressure;  and  3)  there  are  no  heat  and  mass  transfer  factors  resulting  in  temperature  differences  between 
the  liquid  in  the  tank  and  the  vapor  in  the  head  space.  The  behavior  of  the  contents  of  the  tank  while  it 
empties  is  thus  treated  strictly  as  an  equilibrium  process. 

Observations  by  emergency  responders  from  chlorine  railcar  incidents  confirm  that  the  contents  of  the 
tank  auto-refrigerate  leaving  liquid  and  vapor  inside  the  tank  at  ambient  pressure  of  101.325  kPa  (1  atm) 
and  the  liquid  boiling  temperature.  The  final  tank  conditions  are  comprised  of  the  remaining  liquid  mass 
plus  the  vapor  in  the  head  space,  which  is  in  equilibrium  at  the  boiling  temperature.  This  paper  does  not 
address  processes  within  the  tank  subsequent  to  release  of  the  tank  contents  and  auto-refrigeration. 

It  is  assumed  that  there  is  always  liquid  inside  the  hole  until  the  vapor-liquid  interface  reaches  the  top  of 
the  hole.  Under  the  attack  conditions  considered,  it  is  expected  that  there  is  insufficient  energy  imparted 
to  the  tank  contents  to  lead  to  foaming.  The  liquid  head  pressure  is  expected  to  prevent  bubble  formation 
at  the  hole  position  until  the  vapor-liquid  interface  gets  very  close  to  the  top  of  the  hole.4  When  the  vapor 
phase  is  finally  exposed  to  the  atmosphere,  the  head  space  rapidly  depressurizes  from  the  equilibrium 
pressure  to  ambient  pressure.  It  is  during  this  depressurization  that  the  auto-refrigeration  occurs.  Vapor  is 
generated  through  flashing  until  the  liquid  and  vapor  in  the  tank  cool  to  the  boiling  temperature.  Unless 
the  hole  is  a  very  long  vertical  gash,  the  time  between  the  interface  moving  from  the  top  of  the  hole  to  the 
bottom  should  be  small  compared  to  the  overall  tank  empty  time,  so  this  intermediate  phase  is  not 
considered. 
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Tank  volume,  Vtanh  is  estimated  from  initial  mass,  Mtank ,  and  liquid  density,  p^ank'- 


V tank  Mtank/  Pl  ,tank  [l] 

This  volume  is  equal  to  the  partial  volumes  of  remaining  liquid  and  head  space  vapor: 

V lank  My/  pi  +  M ), s/ p  y  [2] 

where  pv  and  pL  represents  the  vapor  and  liquid  densities  respectively  at  the  boiling  temperature.  Mr  and 
Mhs  equal  the  recovered  liquid  mass  and  the  head  space  vapor  mass.  Total  toxic  industrial  chemical  mass 
released,  MTjc,  equals  the  combined  liquid  and  vapor  mass  expelled  from  the  tank: 

Mtic  =  M,ank  -  Mr  -  Mhs  [3  ] 

Until  the  liquid  empties  to  the  top  of  the  hole,  the  pressure  remains  at  saturation  for  the  tank  temperature. 
As  soon  as  the  vapor  space  reaches  the  hole  area,  rapid  depressurization  with  the  associated  auto¬ 
refrigeration  begins.  The  tank  conditions  can  be  determined  at  this  point  just  before  depressurization 
begins  and  combined  with  the  initial  tank  conditions  to  estimate  the  jet  velocity  and  flow  rate  described  in 
Section  2.3. 

Just  before  depressurization,  the  sum  of  the  partial  volumes  of  the  liquid  mass,  MLd,  and  vapor  mass,  Mvd, 
still  equals  the  tank  volume: 


V tank  —  Mid/  p  u  +  Myd/  pyd  [4] 

The  depressurization  point  liquid  density,  pLd,  and  vapor  density,  pvd,  correspond  to  the  saturation 
temperature  at  the  reduced  tank  pressure.  The  energy  required  to  create  the  vapor  mass  comes  from 
cooling  the  current  liquid  and  vapor  contents  in  the  tank.  Assuming  the  liquid  heat  capacity,  Cpljic,  and 
vapor  heat  capacity,  Cpv,tic,  are  constant  over  the  temperature  range  (which  they  essentially  are),  the 
energy  balance  can  be  based  on  the  average  liquid  and  vapor  masses.  The  temperature  decrease, 

AT  =  Tair-  Ttank,  is  written  as: 

AT  =  Mvd  (AHvMnk  +  AHvd)/[(Mtank  +  MLd)Cpuic  +  Myd  Cpv,nc]  [5] 

where  Tmr  is  the  ambient  temperature,  Ttank  is  the  tank  temperature,  AHv  umk  is  the  heat  of  vaporization  at 
the  tank  temperature,  and  AHvd  is  the  heat  of  vaporization  at  the  tank  temperature  prior  to 
depressurization.  Since  additional  mass  vaporizes  as  the  tank  depressurizes,  the  masses  of  vapor  and 
liquid  cannot  be  determined  from  just  Equations  [4]  and  [5]. 

The  depressurization  occurs  adiabatically,  so  enthalpy  is  conserved,  and  the  energy  and  mass  balances 
are: 


MLd  HLd  +  MVd  H Vd  -  MR  HL  +  MVc  (HVd  +  Hv)/2  +  Mhs  Hv  [6] 

Mid  +  Myd  =  Mr+  Myc  +  Mhs  [7] 

H j  d  and  Hvd  are  the  enthalpies  of  liquid  and  vapor  at  the  depressurization  temperature,  and  HL  and  Hv  are 
the  enthalpies  of  liquid  and  vapor  at  the  boiling  temperature.  These  equations  include  both  the  vapor 
leaving  the  tank  as  part  of  the  source  cloud,  MVc,  and  the  vapor  remaining  in  the  head  space  in  order  to 
represent  the  closed  system  of  vapor  and  liquid  masses.  The  temperature  of  the  vapor  leaving  the  tank 
decreases  from  the  tank  temperature  at  the  point  of  depressurization  to  the  boiling  temperature  when  auto- 
refrigeration  is  complete.  For  constant  vapor  heat  capacity,  the  average  of  the  vapor  enthalpies  is  used. 
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The  PV  work  terms  relevant  to  Equations  [5]  and  [6]  are  negligible  compared  to  the  enthalpy  terms.  This 
set  of  equations  must  be  solved  iteratively.  The  procedure  is: 

(1)  Select  AT  and  determine  pLd,  p  vd,  HLd,  and  HVd  from  saturation  tables. 

(2)  Rearrange  Equation  [4]  to  solve  for  MLd,  substitute  it  into  Equation  [5],  and  solve  for  the 
vapor  mass  before  depressurization: 

Myd  ( H tank  V tank  PLd)  CpLJlC  /[(AHV  tank  ~^~AHvd)/  AT  ~^~CpL,TIC  P  Ld  /  PVd  E 'pV.TIC ]  [^] 

(3)  Solve  Equation  [4]  using  MVd  for  the  liquid  mass  in  the  tank  before  depressurization. 

(4)  Solve  Equation  [7]  for  the  vapor  mass  leaving  the  tank. 

(5)  Compute  the  total  enthalpy  before  depressurization,  Hd\ 

Hd  =  MLdHLd  +  MVdHvd  [9] 

(6)  Compute  the  total  enthalpy  after  auto-refrigeration,  Htank\ 

Htcmk  =  MrHl+  MVc  (Hrd  +  Hv)/2  +  Mhs  Hv  [10] 

(7)  If  Hd  is  less  than  H„mk  lower  the  temperature  and  repeat,  and  vice  versa,  until  converged. 

2.3  Velocity  and  Flow  Rate 

Experiments  with  releases  of  superheated  liquids  through  nozzles  have  shown  that  even  well  into  the 
superheated  range  (up  to  around  20  °C),  large  droplets  (several  hundred  microns  in  diameter)  form 
through  breakup  mechanisms.5'9  Ammonia  droplets  larger  than  200  microns  in  diameter  have  been 
calculated  to  be  large  enough  to  reach  the  ground  without  evaporating  from  release  heights  of  several 
meters.5  Since  chlorine  has  the  same  boiling  temperature,  droplet  evaporation  rates  are  comparable.  At 
very  high  superheat  (greater  than  40  °C)  small  droplets  form.  These  droplets  tend  to  be  on  the  order  of 
20-30  microns  in  diameter.5, 9  Rainout  of  large  droplets  at  high  superheat  has  been  found  to  be  minimal, 
suggesting  that  the  fraction  of  larger  droplets  is  very  small.  The  transition  from  formation  of  large 
droplets  with  significant  rainout  to  formation  of  small  droplets  with  minimal  rainout  occurs  around  20  °C 
to  40  °C  superheat.  It  is  assumed  that  the  derivations  here  are  applicable  down  to  20  °C  of  superheat,  with 
the  realization  that  total  cloud  mass  decreases  due  to  rainout  of  larger  droplets. 

For  typical  ambient  wind  conditions,  droplets  smaller  than  10-20  microns  do  not  deposit  onto  surfaces 
due  to  gravitational  settling  because  small  turbulence  eddies  keep  them  airborne.  Droplets  of  all  sizes  still 
deposit  due  to  inertial  and  diffusion  forces.  It  is  likely  that  the  liquid  aerosol  will  coat  the  ground  surface 
and  vegetation  under  the  cloud,  but  that  much  of  the  aerosol  will  not  settle  before  it  evaporates  from  the 
energy  added  by  air  entrained  at  the  cloud  top  surface.  It  is  assumed  that  agglomeration  of  aerosol  to  form 
larger  droplets  is  not  significant.  Any  liquid  coating  the  surface  will  evaporate  quickly  once  the  cloud 
from  which  it  came  is  gone. 

At  storage  pressures  of  5-1 1  atm,  the  experimentally  measured  velocity  of  the  jet  exiting  a  nozzle,  v,  is 
approximately  20  m/s.9  The  relevant  equation  for  the  liquid  flow  is  the  Bernoulli  equation  which  relates 
jet  kinetic  energy  to  tank  pressure,  Ptank,  plus  liquid  head  pressure,  Pi,ead  =  PtankgAh ,  minus  ambient 
pressure,  Patm : 
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Ptank  P head  P atm  P  tank  V  A (2  C-  j)  ) 


[11] 


where  g  is  the  gravitational  acceleration  constant,  and  Ah  is  the  height  between  the  hole  and  the  vapor- 
liquid  interface.  The  orifice  discharge  coefficient,  CD,  accounts  for  frictional  losses  within  the  hole  and  is 
normally  set  to  0.6.  The  tank  fluid  density,  ptank,  is  expected  to  be  that  of  the  liquid  at  the  tank 
temperature;  however,  if  vapor  nucleation  occurs  prior  to  the  orifice,  the  effective  two-phase  density  can 
be  used.  For  the  incidents  considered  here,  the  tank  pressure  and  liquid  density  both  decrease  between  the 
time  of  hole  formation  and  the  time  at  which  the  vapor-liquid  interface  reaches  the  hole.  The  tank 
equilibrium  calculations  in  Section  2.2  predict  that  these  values  do  not  change  greatly. 

Jet  velocity  according  to  the  Bernoulli  equation  is  independent  of  hole  area,  as  long  as  the  discharge 
coefficient  is  not  altered.  The  Bernoulli  equation  assumes  that  the  thickness  of  the  hole  is  negligible.  A 
chlorine  railcar  tank  is  %  inch  thick  and  has  4  inches  of  insulation  and  a  %  inch  thick  metal  jacket;  the 
hole  length  is  therefore  4.875  inches  (0.124  m).  A  hole  in  a  chlorine  railcar  thus  has  a  significant 
thickness,  so  the  jet  experiences  friction  losses  between  the  opening  inside  the  tank  and  the  opening  to 
ambient  air.  The  degree  to  which  the  hole  is  jagged  or  smooth  can  vary  significantly  for  different 
breaching  scenarios,  as  can  the  uniformity  of  the  hole  diameter  across  the  three  layers  in  a  railcar.  For  the 
purposes  of  the  calculations  in  this  work,  the  hole  interior  surface  is  assumed  to  be  smooth  with  a  constant 
diameter.  In  this  case  friction  losses  for  flow  through  a  pipe  apply.  The  friction  factor,/,  is  defined  as: 

/=  AP  d/(2  ptank  v2  L)  [12] 

where  d  is  the  hole  diameter,  L  is  the  hole  length,  and  A P  represents  the  loss  in  pressure  along  the  length 
of  the  hole. 

Flow  through  a  pipe  is  considered  fully  turbulent  for  Reynolds  numbers,  Re  =  d  v  pta„k/pTic,  greater  than 
10,000,  where  jUnc  is  the  viscosity  of  the  liquid  in  the  tank.  The  friction  factor  for  smooth  pipes  is  given 
by  the  Blasius  equation, /=  0.079  Re'0'25,  while  the  friction  factor  for  the  roughness  associated  with 
commercial  pipe  is  slightly  different:10 


/=  0.04  Re016  [13] 

Equating  Equation  [12]  with  Equation  [13]  and  rearranging  gives  the  pressure  drop  along  the  length  of  the 
hole: 


AP  =  0.08  L  ptaJ'84  v1'84  p tic  '16 Ad1 16  [14] 

This  pressure  drop  is  subtracted  from  the  tank  pressure  in  the  Bernoulli  equation  to  give  the  jet  velocity 
adjusted  for  frictional  losses  through  the  length  of  the  hole: 

V  =  CD  [2  (P tcmk  +  P head  —AP  -  P 'at„j/ Ptank]0  5  P  5] 

These  equations  must  be  solved  iteratively.  Note  that,  if  the  pressure  in  the  tank  at  the  hole  equals  the 
chemical  saturation  pressure,  even  this  small  pressure  drop  leads  to  flashing  before  the  orifice.  In  this 
case,  Equations  [12]  through  [15]  do  not  apply;  only  Equation  [1 1]  is  used  with  the  appropriate  fluid 
density.  What  leads  to  frictional  losses  in  the  pipe  for  liquid  flow  now  determines  the  fraction  of  liquid 
flashing  before  the  orifice. 

The  flow  rate,  Q,  for  the  jet  becomes  a  simple  function  of  hole  area: 


Q=  Prank  v  nd2/4 


[16] 
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Tank  empty  time,  tj,  is: 


tf—M tic/ Q  [17] 

Vapor  mass  from  depressurization  and  auto-refrigeration  has  been  included  in  this  calculation,  as  it 
comprises  only  a  small  fraction  of  the  total  chemical  mass  released. 

2.4  Chlorine  Expansion  to  Ambient 

As  the  liquid  passes  through  the  hole,  some  liquid  vaporizes  to  equilibrate  at  ambient  pressure.  The 
remaining  liquid  breaks  up  into  aerosol  and  droplets  due  to  either  mechanical  or  vapor  nucleation  forces. 
The  vapor  expansion  is  adiabatic,  as  no  energy  is  added  to  the  chlorine  as  it  passes  through  the  hole.  No 
work  is  done  on  the  liquid,  either,  so  enthalpy  is  conserved.  The  vapor  expansion  outside  of  the  hole 
occurs  in  all  directions,  although  some  fluid  acceleration  in  the  jet  direction  may  occur  due  to 
backpressure  from  within  the  orifice.  As  long  as  liquid  remains  during  the  adiabatic  expansion,  the  cloud 
forms  at  saturation  temperature  at  the  ambient  pressure  (i.e.,  the  boiling  temperature).  These  conditions 
apply  just  outside  of  the  hole  and  prior  to  the  ambient  air  entrainment  addressed  in  Section  2.5. 

The  fraction  of  vapor  in  the  two-phase  jet  can  be  determined  for  an  arbitrary  starting  mass  of  chemical:6 

MTic  =  Ml+Mv  [18] 

M tic:  Hnc  =  Aft  Hl  +  MyHy  [19] 

Hnc  represents  the  enthalpy  at  the  tank  temperature.  ML  and  My  represent  the  resulting  liquid  and  vapor 
masses  following  adiabatic  expansion.  The  liquid  and  vapor  enthalpy  values  are  at  the  boiling  temperature 
as  used  in  Section  2.2.  The  chemical  mass  represents  a  control  volume  mass  passing  from  inside  the  tank 
to  outside  the  tank.  Note  that  this  derivation  maintains  that  the  fluid  velocity  does  not  change.  Defining  a 
vapor  mass  fraction  in  the  jet,  Fj  =  Mv/MTic,  Equations  [18]  and  [19]  can  be  combined: 


FJ=(Hnc-HI)/(Hv-HI) 


[20] 


Contributions  from  py  volumetric  expansion  and  the  surface  tension  energy  of  the  small  droplets  are 
negligible. 

The  volume  of  mass  released,  Vnc,  is  the  sum  of  the  liquid  and  vapor  partial  volumes: 

Vtic  =  Ml  / Pl  +  My/ py  [21] 

The  density  of  the  two-phase  jet,  ph  is: 

Pj  =  PlPv/[(1  ~Fj)  pv  +  FjpLJ  [22] 

2.5  Air  Entrainment  into  the  Jet 

After  expanding  to  atmospheric  pressure,  the  jet  interacts  with  the  ambient  air  surrounding  the  hole.  The 
entrainment  ratio,  E  =  Vair/Vric,  is  defined  as  the  volume  of  ambient  air,  Vair,  added  to  the  jet  over  the 
volume  of  the  chemical  vapor  and  liquid.  For  this  application  only  air  mixing  into  the  jet  close  to  the  hole 
is  considered.  The  best  example  involves  an  instantaneous  release  of  the  tank  contents  as  a  large  cloud, 
where  the  air  within  that  cloud  volume  mixes  with  the  chemical  instead  of  being  displaced.  Subsequent 
entrainment  of  air  at  the  cloud  or  jet  edges  as  it  moves  downwind  is  not  relevant  to  the  processes 
addressed  in  this  paper. 
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Ambient  air  contains  water  vapor  in  addition  to  the  nitrogen,  oxygen,  and  other  gases  in  dry  air.  The  mass 
fraction  of  water  vapor  in  ambient  air,  FH2o<  is: 

Fmo  =  RHPS  ,H20  Wm(;/[RH  P s,H20  W H20  +  (P atm~  RH  P s,H2o)  ^ air]  [23] 

where  RH  is  the  relative  humidity,  Ps.mo  is  the  vapor  pressure  of  water,  Wmo  is  the  molecular  weight  of 
water,  and  Wair  is  the  molecular  weight  of  dry  air.  The  partial  densities  of  dry  air,  pair.dn>,  and  water  vapor, 
PH20,  are  determined  using  the  ideal  gas  law: 

P  air, dry  (R atm  RH  Ps.H2o)  WJ[(1 000  g/kg)  R  Tair]  [24] 

PH20  =  RHPs,H2o  WH2O/[(1000  g/kg)  R  Tair]  [25] 

where  R  is  the  ideal  constant.  The  moist  air  density,  pair,  is  the  sum  of  these  partial  densities.  Mixing  is 
assumed  to  be  rapid  enough  and  the  chemical  droplets  small  enough  that  they  do  not  rain  out  of  the  cloud 
before  equilibrium  is  achieved.  For  expected  entrainment  ratios,  the  energy  available  from  the  air  and 
water  is  not  sufficient  to  vaporize  all  the  droplets.  Energy  from  the  air  and  water  is  used  to  convert 
chemical  liquid  into  vapor,  so  the  air  and  water  will  cool.  For  the  small  chemical  droplets  present,  heat 
and  mass  transfer  are  very  rapid,  so  the  droplet  temperature  is  not  depressed  during  this  period  of 
evaporation.  The  mixed  cloud  reaches  temperature  equilibrium  when  all  components  are  at  the  chemical 
boiling  temperature  and  ambient  pressure.  Flowever,  since  the  chemical  vapor  can  no  longer  be  at  the 
equilibrium  vapor  pressure  of  1  atm  (101.325  kPa)  with  entrained  air  present,  vapor-liquid  equilibrium  no 
longer  exists.  The  partial  pressure  of  chemical  vapor  cannot  be  raised,  so  the  only  means  to  achieve 
vapor-liquid  equilibrium  is  to  lower  the  equilibrium  pressure  by  cooling  the  cloud  until  the  chlorine  vapor 
pressure  equals  the  chlorine  vapor  partial  pressure.  The  energy  to  cool  the  cloud  comes  from  evaporating 
yet  more  chlorine  liquid,  using  the  chemical  heat  of  vaporization  to  generate  the  temperature  decrease. 
Since  all  partial  pressures  decrease  with  decreasing  temperature,  additional  chemical  vapor  is  generated  in 
this  closed  system  without  exceeding  ambient  pressure.  This  process  is  called  homogeneous  equilibrium 
and  results  in  a  cloud  temperature  below  the  chemical  boiling  temperature.11 

The  equilibrium  processes  affecting  each  cloud  component  are  somewhat  different.  The  entrained  dry  air 
cools  from  ambient  temperature  to  the  final  cloud  temperature  without  any  change  of  phase.  The  water 
vapor  mass  in  air  cools  from  ambient  temperature  to  the  final  cloud  temperature.  Almost  all  of  the  water 
mass  forms  ice  crystals.  Remaining  water  vapor  simply  cools  from  ambient  temperature  to  the  final  cloud 
temperature.  The  mass  fraction  of  water  vapor  remaining,  Fw,„  is  a  function  of  ice  vapor  pressure,  Psjce,  at 
the  final  cloud  temperature: 


Fwv  ~  P s.icJRH P s,H20  [26] 

Chemical  liquid  and  vapor  in  the  cloud  cools  from  the  boiling  temperature  to  the  final  cloud  temperature. 
The  energy  change,  AU,  for  each  component  becomes: 


A  TJ air.dty  Ff air, dry  Fpyiir  (Tair  Tc) 

AUtf20  —  Fwv  MH20  Cpv,H20  (Tair—  Tc) 

HU ice  =  (1  —  Fwv)  Mh20  /A Hv  h20  +  CpL,  H20  (Tair  ~  TfH2o)  + 
AHfM20  +  Cps,H20  (TfH20  ~  Tc) ] 


[27] 

[28] 
[29] 


At/i  —  ((1  -  Fj)  MTic  —  Mve)  CpL  tic  (Tb  —  Tj 


[30] 
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AUv  —  (FjMtic  +  Mys)  Cpvjic  (Tb  -  Tc) 


[31] 


where  Tc  is  the  final  cloud  temperature,  Tb  is  the  boiling  temperature,  Mair  dry  and  MH2o  are  the  entrained 
dry  air  and  water  masses,  Cp_air  and  Cpv,h2o  are  the  heat  capacities  of  dry  air  and  water  vapor,  AHvH2o  and 
AHf  H2o  are  the  heats  of  vaporization  of  water  and  fusion  of  ice,  Cpl,  H20  and  Cps,H2o  are  the  heat  capacities 
of  water  and  ice,  Tfmo  is  the  freezing  temperature  of  water,  and  MVE  is  the  additional  mass  of  chemical 
vaporized  to  reach  equilibrium. 

The  entrained  dry  air  mass  and  water  mass  can  be  related  to  chemical  cloud  mass  using  volume  as  mass 
divided  by  density: 


air, dry  (1  F H2())  Pair  F  M  Tf(y/  p  j 


Mh20  —  F H20  Pair  E  MflC / P.J 

All  energy  from  cooling  and  water  phase  change  is  used  to  convert  chemical  liquid  to  vapor: 


[32] 

[33] 


Mve  AHV  TIC  -  A  U air f dry  +  A  U H20  +  A  Uice  +  A  Ul  +  A  Uv  [34] 

where  AHvjrc  is  the  heat  of  vaporization  of  the  chemical  at  its  boiling  temperature.  Equations  [27]  through 
[33]  can  be  substituted  into  Equation  [34]  to  put  all  energy  components  in  terms  of  cloud  mass.  The 
equation  can  then  be  solved  for  a  defined  mass  evaporation  ratio,  a  =  Mve/MTic- 

.  a  =  ({(1  -  FH2o)  Cpair  (T air  —  T c)  +  Fm,  Fh20  CpyfH20  (T air~  T c)  +  [35] 

(1  -  Fm,)  F H20  [AHV:H20  +  Cpl ,  a 20  (Tair  —  TfH20 )  +AHfH20  + 

Cps.H2o  (TfH2o~  TC)]}E pai,J Pj+[(1  -  F j)Cpijic  + 

Fj  Cpvjic] (Tb  -  Tc))/[AF[v  tic  +(Cpljic  -  Cpv,nc)(Tb  —  Tc) ] 

The  only  unknown  terms  in  Equation  [35]  are  the  final  cloud  temperature  and  water  vapor  mass  fraction, 
with  the  water  vapor  mass  fraction  determined  at  the  final  cloud  temperature  using  Equation  [26]. 

The  chlorine  masses  can  be  written  as: 


ML  —  (1  —  Fj—  a)  Myic 

[36] 

My  =  (Fj  +  a)  Mfic 

[37] 

The  maximum  entrainment  ratio  occurs  when  the  energy  balance  results  in  all  chemical  liquid 
evaporation.  According  to  these  equation,  this  situation  occurs  when  a  =  1-  Fj.  The  maximum 
entrainment  ratio  is  much  higher  than  assumed  for  the  type  of  incident  considered  here. 

2.6  Cloud  Equilibrium 

Determining  the  final  cloud  temperature  requires  a  separate  equation  for  chemical  vapor  partial  pressure, 
Py. 

Pv  =  (1 000  g/kg)  My  R  Tc/Wnc  vc  [38] 

where  Wnc  is  the  toxic  industrial  chemical  molecular  weight,  and  Vc  is  the  cloud  volume.  At  their  boiling 
temperatures,  the  chlorine  vapor  pressure  from  the  ideal  gas  equation  is  two  percent  less  than  that  from 
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the  saturation  table12  and  the  ammonia  vapor  pressure  is  5  percent  less  than  from  the  saturation  table.10 
Therefore,  for  cloud  temperatures  below  the  chemical  boiling  temperature,  then,  ideal  gas  behavior  is 
assumed  to  be  valid  for  the  vapor  components  of  the  cloud.  At  the  final  cloud  temperature,  the  chlorine 
partial  pressure  equals  the  chlorine  saturation  pressure,  Py  =  Ps,tic- 

The  cloud  component  volumes,  V,  are  determined  using  mass  divided  by  pure  component  density  using 
the  mass  relations  presented  in  Section  2.5,  with  the  densities  of  chemical  vapor,  dry  air,  and  water  vapor 
are  computed  using  the  ideal  gas  equation: 


V i  —  (1  —  Fj—  a)  Mnc/ Pcl  [39] 

Vv  =  (1 000  g/kg)  R  Tc  (Fj  +  a)  MTIC/Patm  WTic  [40] 

Vair  =  (1 000  g/kg)  RTC  (1  -  FH20)  Pair  E  MT1C/ Pj  Patm  Wair  [4  1  ] 

Vmo  =  (1 000  g/kg)  R  Tc  F H20  Pwv  Pair  E  MjjC/ Pj  P am  W H20  [42] 

Pice  F H20  (1  F \\  v)  Pair  E  M  TI(/ pj  Pice  [43] 


where  pCL  is  the  density  of  liquid  chlorine  at  the  final  cloud  temperature,  and  plce  is  the  density  of  ice. 
Cloud  volume  is  given  by  the  sum  of  Equations  [39]  through  [43].  Substituting  this  sum  plus  the  vapor 
mass  represented  by  Equation  [37]  into  Equation  [38]  gives  the  chemical  vapor  partial  pressure  in  terms 
of  final  cloud  temperature  and  mass  evaporation  ratio: 

Py  =  (1 000  R  TJWtiq )  (Fj  +  a)/((l  -Fj-  a)/ pcL  +  1000  R  Tc  (Fj  +  o)/Patm  Wnc  +  [44] 

(Epair  /pj){(1000  R  Tc/Patm)[ ( 1  -  FH20)/Wair  +FH20  F,JWH2o]  +  FH20  (1  -  Fm)/pice}) 

As  with  the  tank  equilibrium  calculations  in  Section  2.2,  there  is  no  direct  solution  for  the  final  cloud 
temperature.  The  iterative  procedure  is: 

(1)  Assume  a  final  cloud  temperature  below  the  chlorine  boiling  point  and  determine  pcj,  Ps.ice, 
and  Ps  tic  from  saturation  tables. 

(2)  Calculate  the  fraction  of  water  vapor  remaining  using  Equation  [26]. 

(3)  Calculate  the  mass  evaporation  ratio  using  Equation  [35]. 

(4)  Calculate  the  chemical  vapor  partial  pressure  using  Equation  [44]. 

(5)  If  the  saturation  pressure  exceeds  the  chemical  vapor  partial  pressure,  decrease  the  final  cloud 
temperature,  and  vice  versa,  and  repeat  until  converged. 

Cloud  mass,  Mc,  equals: 


Mc  =  MT1C  (1  +  E  Pair/ Pj)  [45] 

Cloud  density,  pc,  is  given  by  the  cloud  mass  divided  by  the  cloud  volume: 

Pc  =  (1  +  E  pj  pj)/((l  -  Fj  -  a)/ pcL  +  1000  R  Tc  (Fj  +  a)/Patm  W„c  +  [46] 

(E  Pair/ pj){( 1000  R  Tc/Patm)[ ( 1  -  FH20)/Wair  +  FH2„  Fwx/WH2oJ  +  F H20  (1  -  Fw)/picJ) 

The  corresponding  cloud  concentration,  C,  is  given  by  the  total  chemical  mass  divided  by  the  cloud 
volume: 
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[47] 


C  =  1/((1  -  Fj  -  a)/ pcL  +  1000  R  Tc  (Fj  +  a)/Patm  Wnc  + 
(E  paJ  Pj){(1000  R  Tc/Patm)[ (1  -FH20)/Wair  + 
F H20  Fwv/WH20]  +  Fh20  (1  -  Fm,)/ Pice}) 


2.7  Mist  Pool  Theory 

Modern  dense  gas  hazard  assessment  models  use  the  calculations  in  the  preceding  sections,  with  some 
variation.  Results  applying  these  equations  for  a  recent  chlorine  railcar  incident  at  Macdona,  Texas13  for  a 
corresponding  Transportation  Security  Administration  study  are  comparable  to  those  generated  using  two 
proprietary  dense  gas  models.1  The  above  results  are  included  to  support  the  mist  pool  theory  presented  in 
this  section.  Given  release  of  a  toxic  industrial  chemical  through  a  hole  or  rupture  as  a  two-phase  jet,  all 
existing  dense  gas  models  assume  that  the  two  phase  jet  continues  to  exchanges  mass,  momentum,  and 
energy  with  the  surrounding  ambient  air.  This  continued  air  entrainment  quickly  dilutes  the  cloud  and  its 
density,  so  the  cloud  mass  is  carried  downwind  at  a  rate  determined  by  the  wind  velocity  and  loss  of 
initial  momentum.  If  the  slow  vapor  flux  from  boiling  liquid  remaining  in  the  tank  or  evaporating  from 
pooled  liquid  rained  out  from  the  jet  is  ignored,  the  models  predict  that  as  soon  as  the  tank  empty  process 
is  complete,  the  chemical  quickly  leaves  the  incident  area. 

These  model  predictions  are  in  stark  contrast  to  some  observations  where  chlorine  vapor  has  remained  in 
terrain  depressions,  been  blocked  by  tree  lines,  or  collected  in  valleys  for  extended  periods  of  time.13"17 
These  observations  suggest  that  there  are  conditions  for  which  the  passing  ambient  air  cannot  effectively 
remove  the  chemical  as  it  exits  the  hole.  The  result  is  an  accumulation  of  chemical  vapor  and  aerosol  that 
persists  at  the  incident  location  for  a  measurable  period  of  time  beyond  the  tank  empty  time.  This  cloud  of 
highly  dense  chemical  vapor  and  aerosol  is  defined  as  a  mist  pool.  This  paper  derives  a  criterion  that  can 
be  used  to  estimate  when  a  mist  pool  will  form  and  provides  the  equations  characterizing  the  behavior  of 
a  mist  pool  once  it  has  formed.  The  calculations  are  then  applied  to  the  example  chlorine  and  ammonia 
incidents. 

A  mist  pool  is  proposed  to  form  when  the  rate  of  toxic  industrial  chemical  passing  through  the  hole 
exceeds  the  capacity  of  the  air  to  remove  it  from  the  incident  location.  The  chemical  vapor  and  aerosol 
that  is  not  taken  up  by  passing  air  collects  at  the  incident  site,  continuing  to  add  to  the  accumulated  cloud 
mass  until  the  tank  empty  time  is  reached.  The  high  cloud  density  leads  to  the  slumping  behavior  that 
dense  gas  models  already  treat  for  instantaneous  releases.  The  cloud  radius  expands  over  flat  and 
unobstructed  terrain,  with  considerable  turbulence  within.  Observations  from  field  trials  suggest  the 
visible  cloud  height  stabilizes  at  approximately  1  m. 18-25  The  final  cloud  form  becomes  a  flat  cylinder  of 
large  diameter.  Air  is  entrained  by  the  combined  events  of  interaction  of  the  jet  with  passing  air  and 
turbulent  entrainment  at  the  cloud  edges.  Since  the  entrainment  ratio  with  passing  air  is  reduced,  the 
entrainment  at  the  cloud  edges  is  assumed  to  balance  that  reduction  to  lead  to  an  entrainment  ratio 
comparable  to  the  jet  calculations  in  Section  2.6,  which  is  assumed  to  Equation  [5],  The  cloud,  which 
contains  the  combined  mass  of  chemical  vapor  and  aerosol,  air,  water  vapor,  and  ice,  can  be  represented 
by  a  1  m  thick  disk  with  the  components  uniformly  distributed  throughout: 


Mc  =  npc  Dc  Hc/4 


[48] 


where  Dc  is  the  cloud  diameter,  and  H,  is  the  cloud  height  or  thickness.  If  the  average  mass  flow  rate 
through  the  hole  is  considered  and  it  is  assumed  that  the  cloud  height  is  constant,  the  cloud  diameter  can 
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be  represented  as  a  function  of  time: 

Dc  =  [4  Q  (1  +  E paJ pj)/ n  pc  HJ05  t05  [49] 

where  t  represents  time  from  hole  formation.  Although  the  cloud  diameter  increases  with  t0'5,  the  cloud 
top  area  increases  linearly  with  time.  The  average  cloud  top  area  during  cloud  formation  is  then  just  half 
the  final  cloud  top  area. 

If  a  mist  pool  forms  according  to  the  criterion  proposed  below  in  Section  2.8,  the  cloud  is  not  moved  by 
the  wind.  Rather,  it  remains  stationary  and  provides  an  area  vapor  flux  to  the  passing  air  until  it  is 
depleted.  The  concept  of  a  stationary  area  vapor  source  is  not  new.  The  HGSYSTEM26  and  DEGAD1S27’ 
28  models  address  the  situation  where  the  rate  of  vapor  formation  from  an  evaporating  liquid  pool  exceeds 
the  capability  of  the  passing  air  to  remove  it.  An  expanding  secondary  source  blanket  is  then  formed 
above  the  liquid  pool  and  is  treated  as  a  stationary  vapor  source,  with  vapor  exchange  to  the  passing  air  at 
the  cloud  top.  This  situation  persists  until  vapor  generation  conditions  change.  The  area  vapor  source 
approach  has  not  yet  been  applied  to  the  vapor,  aerosol,  and  entrained  air  mix  from  a  two-phase  jet.  The 
mist  pool  theory  in  this  paper  extends  the  secondary  source  blanket  approach  used  by  HGSYSTEM  and 
DEGAD1S. 

2.8  Mist  Pool  Formation  Criterion 

The  HGSYSTEM  and  DEGADIS  models  both  compute  a  maximum  entrainment  or  take-up  rate  based  on 
the  area  of  the  top  surface  of  the  secondary  source  blanket.  As  long  as  the  mass  flux  from  the  evaporating 
pool  exceeds  the  maximum  entrainment  rate,  the  resulting  cloud  remains  centered  at  the  release  location 
and  increases  in  mass  equal  to  the  difference  between  the  rates.  Once  no  additional  mass  is  being 
generated  or  the  mass  rate  becomes  less  than  the  maximum,  the  models  decrement  the  cloud  mass  by  the 
maximum  entrainment  rate  until  no  vapor  remains.  The  mist  pool  process  proposed  here  applies  the  same 
maximum  entrainment  rate  to  the  jet  from  railcar  incidents.  If  the  mass  flow  rate  through  the  hole  does 
not  exceed  the  maximum  entrainment  rate,  no  stationary  cloud  forms.  In  this  case,  the  existing  modeling 
capability  (e.g.,  in  HGSYSTEM,  DEGADIS,  and  other  modem  dense  gas  models)  is  considered 
appropriate. 

Entrainment  rates  are  computed  by  dense  gas  models  to  estimate  the  rate  of  mass  transfer  across  the 
effective  surface  boundary  between  ambient  air  and  the  cloud  or  jet.  The  typical  use  is  to  determine  the 
rate  of  entrainment  of  ambient  air  into  the  cloud  or  jet  mass;  however,  the  entrainment  rates  are  not 
directional,  so  they  can  just  as  readily  be  used  to  estimate  the  entrainment  of  chemical  from  the  jet  into 
ambient  air.  The  HGSYSTEM  and  DEGADIS  models  contain  similar  entrainment  rate  equations  for  both 
directions,  the  use  of  which  becomes  situation  dependent.  The  vertical  entrainment  velocity  is  a  function 
of  the  Richardson  number,  Ri*: 


Ri*  g  ( pc  Pair)  Pair  U*  [50] 

where  u*  is  the  friction  velocity  for  the  ambient  wind  conditions.  By  itself,  the  Richardson  number 
provides  an  estimate  of  the  relative  resistance  of  the  cloud  mass  to  the  passing  air.  Vertical  entrainment 
velocity,  tv*,  for  the  secondary  source  blanket  in  DEGADIS27  28  is  given  as: 

w*  =  0. 735  u *  (1  +  a)/(0. 88  +  0. 046  Ri*104)  [5 1  ] 

In  this  equation,  a  equals  the  exponent  for  the  vertical  wind  profile  power  law  equation.  Pulling  together 
all  contributions  to  the  vertical  entrainment  velocity  from  HGSYSTEM  simplifies  to:26 

W*  =  2.87  u*/Ri*  [52] 
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for  Ri«  >12.  An  alternate  relation  is  used  for  lower  Richardson  numbers: 


w*  =  0. 70  u*/(l  +  0.8  Ri*)0'5  [53] 

A  third  equation  for  vertical  entrainment  velocity  is  given  by  Briggs29: 

w*  =  0.65  u*/(l  +  0.2  Ri*)  [54] 

The  vertical  entrainment  velocity  equations  were  derived  from  data  on  all-vapor  sources  with  Richardson 
numbers  up  to  about  20;29  given  the  lack  of  experimental  data  and  equations  for  two-phase  entrainment 
velocities,  the  relations  described  above  are  used. 

The  vertical  entrainment  predictions  were  compared  for  three  ambient  wind  conditions,  all  for  a  grassy 
surface  (5  cm  roughness  length)  on  open  terrain:  unstable  conditions  with  a  2  m/s  wind  speed  at  2  m 
above  the  ground  surface,  neutral  conditions  with  a  5  m/s  wind  speed,  and  stable  conditions  with  a  2  m/s 
wind  speed.  Three  conditions  for  the  chlorine  and  ammonia  incidents  were  considered:  cloud  formation 
(Ri*  from  24  to  2087),  no  remaining  aerosol  (Ri*  from  6  to  689),  and  a  Richardson  number  of  zero.  The 
cloud  density  after  all  aerosol  has  evaporated,  pcp  is  obtained  by  subtracting  the  liquid  aerosol  mass  and 
partial  volume  from  Equation  46: 

Per  =(Fj+a  +  E pj pj)/(1000  R  Tc  (Fj  +  cc)/Patm  Wnc  +  [55] 

(E  pj  Pj)  {(1000  R  TJPaJ[(l  -Fmo)/Wair  +  FH20FJWH2O]  +  Fmo  (1  -FJ/pice}) 

The  corresponding  cloud  concentration,  Cy,  is  needed  to  estimate  mass  flux: 

Cf  =  (Fj  +  a)/(1000  R  Tc  (Fj  +  cc)/Patm  WT1C  +  [56] 

(E  pairf  p.j)  {( 1 000  R  Tc/PaJ[(l  -  Fmo)/Wair  + 

F H20  Fwv/Wh2o]  +  Fh20  (7  -  Fm,)/ Pice}) 

The  comparison  of  the  velocities  shows  a  typical  range  of  about  a  factor  of  5  across  the  three  model 
approaches.  The  HGSYSTEM  and  Briggs  equations  produce  very  similar  values  for  all  conditions,  while 
the  DEGADIS  equations  always  result  in  the  highest  velocity.  The  E1GSYSTEM  equations  have  been 
selected  for  mist  pool  vapor  flux  computations. 

The  maximum  mass  entrainment  rate,  Qm  at  the  cloud  top  equals  Qv  =  Af,,  where  f,  is  the  vapor  mass  flux 
and  A  is  the  cloud  top  area.  The  area  equals  half  the  final  cloud  top  area.  Vapor  mass  flux  is  given  by  f,  = 
C  w*.  The  maximum  entrainment  rate  equation  for  Ri*  >12  becomes: 

Qv  =  2.87  C  u*  nDc2/8  Ri*  [57] 

The  proposed  criterion  for  mist  pool  formation  is  that  the  ratio  of  the  maximum  entrainment  rate  over  jet 
flow  rate,  QJQ ,  be  less  than  one. 

Additional  consideration  must  be  made  for  the  fact  that  the  vertical  entrainment  velocity  equations  were 
derived  from  data  taken  over  homogeneous,  flat  terrain  and  do  not  account  for  the  orientation  of  the  jet.  A 
jet  orientation  resists  being  taken  up  by  the  passing  air  when  it  is  pointing  into  the  wind,  or  downward. 
When  the  initial  cloud  behavior  is  dominated  by  gravitational  slumping,  terrain  slope  and  depressions 
affect  cloud  stability.  Any  mass  settling  into  terrain  depressions  resists  being  removed  by  the  air.  Finally, 
obstructions  such  as  buildings  and  tree  stands  provide  a  barrier  to  cloud  movement.  If  buildings  or  trees 
obstruct  the  initial  cloud  movement,  it  resists  being  moved  by  the  wind.  The  effects  of  the  jet  aiming 
upwind  or  downward,  terrain  depressions,  and  building  and  tree  obstructions  are  roughly  approximated  by 
decreasing  the  entrainment  velocity  (and  thus  the  entrainment  rate)  by  the  factor  of  5  found  to  represent 
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the  spread  of  entrainment  velocities  under  normal  open  terrain  conditions.  Enhanced  entrainment  from  the 
jet  aiming  with  the  wind  is  likewise  accounted  for  using  a  factor  of  5  increase.  The  first  conditions  are 
referred  to  here  as  sheltered  and  the  latter  as  enhanced.  Equation  [57]  is  rewritten  as  the  mist  pool 
formation  criterion  with  the  entrainment  factor,  B: 

QJQ  =  2.87  BC Pair  u J  ttDc2/[8  Qg(pc-  pair)  HJ  <  1  [58] 

where  B  =  5  for  enhanced  entrainment,  1  for  normal  entrainment,  and  0.2  for  sheltered  entrainment.  Note 
that  the  ratio  increases  strongly  with  increases  in  friction  velocity  and  final  cloud  diameter  and  decreases 
linearly  with  increases  in  jet  mass  flow  rate  and  cloud  height.  The  cloud  concentration  and  density  are 
related;  the  effect  is  that  the  ratio  changes  slowly  with  air  entrainment  ratio  and  jet  density.  Ratios  close  to 
unity  signify  that  if  a  mist  pool  forms,  it  will  only  last  a  short  while  beyond  the  tank  empty  time.  If  a  mist 
pool  does  not  form,  the  cloud  mass  will  not  be  easily  moved  by  the  passing  wind,  and  parts  of  it  will  be 
retained  within  terrain  depressions  and  blocked  by  obstructions.  A  combination  of  bulk  cloud  movement 
at  reduced  speed  with  a  significant  vapor  flux  into  passing  air  is  likely. 

2.9  Mist  Pool  Vapor  Flux 

The  mist  pool  process  proposed  for  massive  releases  of  toxic  industrial  chemical  follows  three  stages.  The 
criterion  in  Section  2.8  determines  if  a  mist  pool  is  expected  to  form  under  the  release  conditions.  If  a 
mist  pool  does  not  form,  the  normal  dense  gas  modeling  of  the  jet  applies. 

2. 9. 1  Tank  Emptying 

During  the  first  stage,  the  cloud  expands  as  a  disk  of  constant  height  with  a  diameter  given  by  Equation 
[49].  This  stage  continues  until  the  tank  empty  time  is  reached.  For  conditions  meeting  the  mist  pool 
criterion,  this  time  is  several  minutes  or  less.  Chemical  vapor  is  entrained  at  the  cloud  top  during  this 
stage.  It  is  assumed  that  all  vapor  lost  is  replaced  by  evaporated  aerosol,  so  the  cloud  density,  chlorine 
concentration,  and  Richardson  number  all  decrease  somewhat.  The  dynamic  conditions  keep  the  cloud 
well  mixed  throughout. 

The  average  entrainment  rate  during  stage  1,  Qvj,  is  given  by  Equation  [57].  In  order  to  account  for  the 
mass  lost  during  this  stage  the  fraction  of  aerosol  mass  remaining,  1  -  Fj-  a  -  Qvi/Q ,  is  determined.  The 
cloud  density,  pcl,  and  concentration,  C/,  are  recalculated  for  this  additional  aerosol  mass  lost: 


Pci  =  ( 1 +  EpJpj  -  QJQ) /((l  -Fj-  a-  QJQ)/ PcL  +  1000  R  Tc  (Fj  +  q)/Patm  W„c  +  [59] 

(E  paJ  Pj)  {(1000  R  TJP atm)[ ( 1  -  FH20)/Wair  +  FH20  FJWmo]  +  FH20  (1  -  FJ/pJ) 

Ci  =(1  -  QJQ)/((1  -Fj- a-  QJQ)/ pcL  +  1000 R  Tc  (Fj  +  a)/Patm  W„c  +  [60] 

(E  paJ  pj)  {(1000  R  TJPatm)[(l  -  FH20)/Wmr  +  Fmo  FJWmo]  +  FH20  (1  -  Fj/pj) 

The  Richardson  number  and  average  entrainment  rate  are  then  recalculated  using  Equations  [50]  and  [57], 
respectively.  The  aerosol  mass  remaining  at  the  end  of  stage  1  is  then  calculated  as  ML1  =  (1  -  Fj-  a  - 
QJQ)  Mnc-  If  the  aerosol  mass  remaining  becomes  less  than  zero  for  either  this  equation  or  the  estimate 
before  recalculation,  the  aerosol  evaporates  as  the  tank  empties,  and  the  average  entrainment  rate  is  reset 
to  Qvi  =  (1  -  Fj- a)  Q.  Chemical  mass  at  the  end  of  stage  1  equals  Mi  =  MT]C  -  Qvi  tf 
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2.9.2  Aerosol  Evaporation 


During  the  second  stage,  the  cloud  maintains  a  constant  diameter  and  generates  vapor  from  the  cloud  top 
surface.  The  vapor  flux  is  replaced  by  evaporating  aerosol.  If  the  aerosol  evaporates  as  the  tank  empties, 
this  stage  is  skipped.  Any  air  mass  lost  with  the  vapor  at  the  top  surface  is  replaced  with  an  equivalent 
amount  of  moist  air.  The  air  substitution  plus  heat  flux  from  the  ground  surface  provide  the  energy  for 
evaporating  the  aerosol,  and  all  components  are  assumed  to  be  well  mixed  at  all  times.  This  assumption 
should  be  valid  for  several  minutes  beyond  the  tank  empty  time,  after  which  the  internal  cloud  turbulence 
subsides,  and  aerosol  begins  to  deposit  onto  the  ground  surface.  Mass  deposited  as  a  film  on  the  ground 
surface  below  the  cloud,  though,  and  then  evaporates  as  part  of  the  all-vapor  cloud  in  the  next  stage.  The 
vapor  flux  is  assumed  constant  during  this  step.  As  long  as  aerosol  remains,  the  cloud  temperature  does 
not  change.  Cloud  volume  decreases  slightly  with  time  due  to  loss  of  the  aerosol,  but  this  insignificant 
change  is  ignored.  The  cloud  density,  chlorine  concentration,  and  Richardson  number  decrease  steadily  as 
aerosol  evaporates. 

The  entrainment  rate  at  the  end  of  stage  1  equals  twice  the  recalculated  value  from  Equation  [57].  The 
entrainment  rate  at  the  end  of  stage  2  is  then  calculated  as: 

Q,f  =  2.87  B  Cfu*  7tDc2/4  Ri*f  [61] 

where  Qvl  is  the  entrainment  rate  when  no  aerosol  remains  and  Rij is  computed  using  the  cloud  density 
after  all  aerosol  has  evaporated.  The  entrainment  rate  during  stage  2,  Qv2,  is  computed  as  the  average  of 
gv/and  that  at  the  end  stage  1.  The  duration  of  stage  2,  t2,  equals  t2  =  ML1/Qv2,  and  the  chemical  mass 
remaining  at  the  end  of  stage  2  becomes  M2  =  (Fj  +  a)  MTic- 

2. 9. 3  Vapor  Depletion 

The  third  stage  begins  when  all  aerosol  mass  has  evaporated,  but  it  cannot  begin  before  the  tank  empty 
time.  The  vapor  mass  is  assumed  to  be  stripped  from  the  cloud  top  with  no  air  entrainment  into  the  cloud. 
The  result  is  a  vapor  cloud  of  constant  composition  that  is  linearly  decreasing  in  thickness.  Cloud  density 
and  concentration  remain  constant,  but  the  Richardson  number  decreases  with  the  cloud  thickness.  The 
starting  vapor  mass  for  stage  3  always  begins  with  Equation  [68].  During  stage  3,  the  Richardson  number 
crosses  a  value  of  12,  so  two  entrainment  equations  are  needed.  The  entrainment  rate  when  the 
Richardson  number  equals  12,  Qv3a,  is: 


Qv3a  =  2.87  B  Cfu*  kDc2/(4)(12)  [62] 

The  fraction  of  the  cloud  height,  fc,  when  the  Richardson  number  equals  12  is  determined  next: 

fc  Pair  U*  (i^)/§(pcf  Pair)  i ! c  [63] 

The  duration  of  the  first  part  of  stage  3  before  the  Richardson  number  reaches  12  is  t3a  =  2  (1  -  f  )  M2/(Qvf 
+  Qv3a)-  The  rate  when  the  last  vapor  mass  is  stripped  off,  Qv3b,  is  determined  for  a  Richardson  number  of 
0: 


Qv3h  =  0.70  BCu*  nD2/4  [64] 

The  duration  of  the  second  part  following  the  Richardson  number  reaching  12  is 

hb  =  2  fc  M2/(Qv3a  +  Qv3b)- 
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2.9.4  Combined  Vapor  Flux 


The  total  mist  pool  duration,  tent,  which  includes  the  tank  empty  time,  equals  the  sum  of  four  times: 


tent  tf  t  2  tja  + 


[65] 


The  average  entrainment  rate  equals  Qent  =  MTIC  /tem.  Regardless  of  the  how  long  the  mist  pool  lasts 
beyond  the  tank  empty  time,  the  source  geometry  changes  dramatically  from  a  narrow  two-phase  jet  to  a 
large  diameter,  flat  cloud.  Note  that  the  entrainment  rate  increases  greatly  as  the  mist  pool  progresses 
from  cloud  formation  to  all  vapor.  Based  on  the  very  high  mass  entrainment  rates  calculated  when  only 
vapor  remains,  the  friction  velocity  will  have  to  be  very  low  in  order  for  a  secondary  vapor  blanket  to 
form.  The  air  will  have  to  be  stable  and  still.  With  addition  of  the  mist  pool  aerosol,  resistance  to 
entrainment  at  the  cloud  top  is  greatly  increased,  so  a  mist  pool  can  form  under  normal  wind  conditions 
for  a  high  jet  flow  rate. 

3.0  RESULTS  AND  DISCUSSION 

Two  example  incidents  are  used  to  demonstrate  the  mist  pool  concept;  both  involve  a  standard  90-ton 
railcar.  The  railcar  volume  is  filled  with  chlorine  and  then  anhydrous  ammonia.  It  is  assumed  that  the 
ambient  temperature  and  tank  temperature  are  20  °C  (293. 15  K).  A  full  chlorine  railcar  normally  has  a 
small  vapor  head  space  or  ullage  of  about  10  percent  of  the  tank  volume.  To  simplify  calculations,  it  is 
assumed  that  there  is  no  vapor  head  space  at  the  time  of  hole  formation.  The  hole  is  assumed  to  be 
positioned  about  1/3  of  the  way  up  from  the  tank  bottom  to  the  top,  so  only  liquid  is  present  on  the  inside 
at  the  hole  during  the  release.  The  final  recovered  chlorine  liquid  mass  is  20  tons;  the  recovered  ammonia 
mass  is  the  equivalent  volume.  Initial  chlorine  mass  is  81,596  kg,  and  recovered  chlorine  liquid  mass  is 
18,132  kg.  The  initial  mass  of  ammonia  is  determined  by  relative  liquid  densities  at  20  °C.  Physical 
properties  for  chlorine  and  ammonia  are  provided  in  Table  1.  The  initial  ammonia  mass  is  35,328  kg.  The 
liquid  densities  at  the  final  tank  temperature  are  used  to  determine  the  recovered  mass  of  ammonia  to  be 
7931  kg.  Tank  volume  (Equation  [1])  is  57.93  m3.  Head  space  vapor  mass  (Equation  [2])  is  171  kg  of 
chlorine  vapor  and  42  kg  of  ammonia  vapor.  The  chlorine  incident  releases  63,292  kg,  and  the  ammonia 
incident  releases  27,355  kg. 

Table  1.  Physical  Properties  of  Ammonia  and  Chlorine 


|  PHYSICAL  PROPERTY 

VARIABLE 

CHLORINE10  1Z’ 30 

AMMONIA10  [ 

boiling  temperature  (C/K) 

Tb 

-34.0/239.15 

-33.4/239.75 

molecular  weight  (g/mol) 

Wnc 

70.91 

17.03 

liquid  density  at  20  °C  (kg/m3) 

PL.tank 

1408.6 

609.9 

liquid  density  at  Tb  (kg/m3) 

Pl 

1557.9 

681.4 

vapor  density  at  Tb  (kg/m3) 

Pv 

3.6985 

0.9075 

heat  of  vaporization  at  20  °C  (kJ/kg) 

AH Vi  tank 

255.13 

1188.26 

liquid  heat  capacity  (kJ/kg-K) 

CpL,TIC 

0.926 

4.65 

vapor  heat  capacity  (kJ/kg-K) 

CpV.TIC 

0.50 

2.07 

liquid  enthalpy  at  Tb  (kJ/kg) 

Hl 

237.64 

29.4 

vapor  enthalpy  at  Tb  (kJ/kg) 

Hv 

525.72 

1401 

liquid  enthalpy  at  20  °C  (kJ/kg) 

Htic 

288.81 

275.19 

liquid  viscosity  (cP) 

Ptic 

* 

0.12 

no  data 


15 


3.1  Jet  Characterization 


The  tank  emptying  equilibrium  results  (Section  2.2)  for  the  chlorine  incident  converge  at  a  tank 
temperature  before  depressurization  of  289.0  K  (15.9  °C),  or  a  decrease  in  temperature  of  4.1  °C.  The 
degree  of  superheat  decreases  from  54  °C  to  50  °C,  so  it  remains  very  high.  The  pressure  in  the  tank 
decreases  from  677.9  kPa  (6.69  atm)  to  606.9  kPa  (5.99  atm),  and  the  vapor  density  decreases  from  21.63 
kg/m3  to  19.49  kg/m3.  The  liquid  density  increases  slightly  from  1408.6  kg/m3  to  1420.8  kg/m3.  The 
results  for  the  ammonia  incident  converge  at  a  tank  temperature  before  depressurization  of  290.4  K  (17.3 
°C),  or  a  decrease  in  temperature  of  2.8  °C.  The  degree  of  superheat  decreases  from  53  °C  to  50  °C.  The 
pressure  in  the  tank  decreases  from  861.3  kPa  (8.50  atm)  to  792.4  kPa  (7.82  atm),  and  the  vapor  density 
decreases  from  6.72  kg/m3  to  6.20  kg/m3.  The  liquid  density  increases  slightly  from  609.9  kg/m3  to  613.9 
kg/m3. 

A  single  round  hole  of  0. 10  m  diameter  (4  inches)  is  considered  for  a  release  from  a  railcar.  The  pressure 
drop  through  the  hole  (Equation  [14])  is  5.5  kPa  (0.054  atm)  for  the  chlorine  incident  and  7.7  kPa  (0.076 
atm)  for  the  ammonia  incident.  The  average  chlorine  jet  velocity  (Equation  [15])  is  18.9  m/s,  and  the 
average  ammonia  jet  velocity  becomes  34.8  m/s.  The  final  chlorine  Reynolds  number  of  1.49E5  and 
ammonia  Reynolds  number  of  1.05E5  are  both  still  well  into  the  turbulent  flow  regime.  Both  jet  velocities 
are  comparable  to  the  value  of  20  m/s  cited  in  Section  2.3.  Average  chlorine  flow  rate  (Equation  [16])  is 
166  kg/s,  and  the  tank  empty  time  (Equation  [17])  is  382  s  (6.4  min).  The  average  ammonia  flow  rate  is 
121  kg/s,  and  the  tank  empty  time  is  226  s  (3.8  min).  For  the  position  of  the  hole  in  the  railcar  considered 
here,  the  initial  liquid  head  pressure  is  about  1 1  to  25  kPa  (0.1 1  to  0.25  atm),  so  the  pressure  at  the  orifice 
does  not  reach  saturation  until  the  vapor-liquid  interface  approaches  the  top  of  the  hole. 

The  chlorine  incident  involves  expulsion  of  59,290  kg  of  liquid  before  depressurization  and  4002  kg  of 
vapor  during  depressurization  and  auto-refrigeration.  The  ammonia  incident  involves  expulsion  of  25,628 
kg  of  liquid  before  depressurization  and  1727  kg  of  vapor  during  depressurization  and  auto-refrigeration. 
Therefore,  94%  of  the  total  chemical  mass  in  each  cloud  is  released  prior  to  depressurization. 

For  the  adiabatic  expansion  of  the  jet  exiting  the  hole  discussed  in  Section  2.4,  both  the  chlorine  and 
ammonia  jets  are  18%  vapor  and  82%  liquid.  The  density  of  the  chlorine  jet  (Equation  [22])  is  20.8 
kg/m3,  and  the  density  of  the  ammonia  jet  is  5.1  kg/m3.  The  adiabatic  expansion  leads  to  an  increase  in 
volume  of  68  times  that  of  the  chlorine  liquid  inside  the  tank  and  121  times  that  of  the  ammonia  liquid 
inside  the  tank. 

For  the  air  entrainment  addressed  by  Section  2.5,  physical  property  values  and  constants  needed  for 
equations  dealing  with  water  and  air  are  included  in  Table  2.  Relative  humidity  is  assumed  to  be  0.50. 

The  resulting  water  mass  fraction  (Equation  [23])  is  0.0072.  At  20  °C  and  50%  relative  humidity,  the 
partial  density  of  dry  air  (Equation  [24])  equals  1.1900  kg/m3  and  the  partial  density  of  water  vapor 
(Equation  [25])  equals  0.0086  kg/m3.  The  moist  air  density  then  equals  1.1986  kg/m3. 

The  final  cloud  properties  described  in  Section  2.6  have  been  determined  for  an  entrainment  ratio  of  5. 
The  final  chlorine  cloud  temperature  is  -57  °C,  and  the  saturation  pressure  is  33.0  kPa  (0.325  atm).  The 
fraction  of  water  vapor  remaining  is  0.0091.  16%  of  the  total  chemical  cloud  mass  evaporates  in  addition 
to  the  initial  18%  from  adiabatic  expansion  of  the  jet.  The  final  ammonia  cloud  temperature  is  -54  °C,  and 
the  saturation  pressure  is  32.0  kPa  (0.3 16  atm).  The  fraction  of  water  vapor  remaining  is  0.01 14,  and  14% 
of  the  cloud  mass  evaporates  in  addition  to  the  initial  18%.  Achieving  both  temperature  and  vapor-liquid 
equilibrium  following  initial  air  entrainment  results  in  both  clouds  cooling  by  more  than  20  °C  below  the 
chemical  boiling  temperatures.  The  chlorine  incident  cloud  density  (Equation  [46])  is  4.9  kg/m3,  or  4.1 
times  that  of  ambient  air.  The  ammonia  incident  cloud  density  is  2.0  kg/m3,  or  1.7  times  that  of  ambient 
air.  The  chlorine  concentration  (Equation  [47])  is  3.8  kg/m3,  and  the  ammonia  cloud  concentration  is  0.9 
kg/m3.  The  chlorine  cloud  density  after  all  aerosol  has  evaporated  (Equation  [55])  is  2.4  kg/m3,  and  the 
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corresponding  ammonia  cloud  density  is  1.4  kg/m3.  The  chlorine  and  ammonia  cloud  concentrations  after 
aerosol  evaporation  (Equation  [56])  are  1.3  kg/m3,  and  0.3  kg/m3,  respectively. 

Table  2.  Other  Physical  Properties  and  Contents 


I  PHYSICAL  PROPERTY/CONSTANT 

VARIABLE 

VALUE10’31  1 

water  vapor  pressure  at  20  C  (kPa) 

P  S.H20 

2.3 

ice  vapor  pressure  at  -10  C  (kPa) 

P  . 

•  s,ice 

0.27 

molecular  weight  of  water  (g/mol) 

1/14/20 

18.016 

molecular  weight  of  dry  air  (g/mol) 

Wair 

28.96 

water  freezing  temperature  (C/K) 

Tf:H20 

0/273.15 

ideal  gas  constant  (m3-kPa/mol-K) 

R 

0.008314 

gravitational  acceleration  constant  (m/s2) 

9 

9.81 

Pi 

P 

3.14159 

heat  capacity  of  air  (kJ/kg-K) 

6'pa//- 

1.05 

heat  capacity  of  ice  (kJ/kg-K) 

CPS,H20 

1.97 

heat  capacity  of  water  (kJ/kg-K) 

CpL,H20 

4.20 

heat  capacity  of  water  vapor  (kJ/kg-K) 

Cpv,H20 

1.88 

water  heat  of  fusion  (kJ/kg) 

DHfH20 

333.7 

water  heat  of  vaporization  (kJ/kg) 

DHvh20 

2261 

density  of  ice  (kg/m3) 

r ice 

915 

3.2  Mist  Pool  Behavior 

For  aim  cloud  height,  the  chlorine  incident  mist  pool  cloud  diameter  as  a  function  of  time  (Equation 
[49])  becomes  Dc  =  7.41  t"'\  with  a  final  diameter  of  145  m.  The  ammonia  cloud  diameter  as  a  function 
of  time  is  Dc  =  12.82  t"  '\  with  a  final  diameter  of  193  m.  The  mist  pool  formation  criterion  derived  in 
Section  2.8  has  been  applied  to  the  chlorine  and  ammonia  clouds.  The  entrainment  ratio  (Equation  [58]) 
for  the  chlorine  incident  under  neutral  conditions  (friction  velocity  of  0.54  m/s)  and  sheltered  entrainment 
is  0.56,  and  the  ratio  for  the  ammonia  incident  is  1.48,  so  a  mist  pool  is  expected  to  form  for  the  chlorine 
release  but  not  for  the  ammonia  release.  The  ratios  for  both  incidents  are  within  a  factor  of  two,  so 
relatively  small  changes  in  the  release  conditions  can  switch  the  expected  behavior. 

Since  the  ammonia  incident  does  not  meet  the  mist  pool  formation  criterion,  the  vapor  flux  calculations  in 
Section  2.9  can  only  be  done  for  the  chlorine  incident.  For  the  tank  emptying  stage  defined  in  Section 
2.9.1,  the  initial  chlorine  entrainment  rate  is  93  kg/s.  The  fraction  of  aerosol  remaining  becomes  0.10.  The 
cloud  density  and  concentration  are  recalculated  to  be  2.8  kg/m3  and  1.7  kg/m3,  respectively.  The 
Richardson  number  decreases  from  105  to  34,  and  the  average  entrainment  rate  is  recalculated  to  be  96 
kg/s.  The  fraction  of  aerosol  remaining  is  updated  to  be  0.08,  and  the  aerosol  mass  remaining  in  the  cloud 
at  the  end  of  stage  1  is  4965  kg.  Chlorine  mass  in  the  cloud  is  26,472  kg.  The  average  entrainment  rate  for 
the  aerosol  evaporation  stage  (Section  2.9.2)  is  194  kg/s,  with  a  duration  of  26  s.  The  chemical  mass 
remaining  is  21,508  kg.  The  entrainment  rate  for  the  first  part  of  the  vapor  depletion  stage  (Section  2.9.3) 
is  376  kg/s  with  a  duration  of  37  s.  The  fraction  of  the  cloud  height  when  the  Richardson  number  reaches 
12  is  0.35,  so  65%  of  the  vapor  mass  is  stripped  off  during  the  first  part  of  that  stage.  The  entrainment  rate 
for  the  second  part  of  the  vapor  depletion  stage  is  1092  kg/s  with  a  duration  of  7  s.  The  total  mist  pool 
duration  for  the  chlorine  incident  is  452  s  (7.5  min),  with  an  average  entrainment  rate  of  140  kg/s.  Since 
the  tank  empty  time  is  382  s,  the  mist  pool  adds  only  70  s  to  the  cloud  duration  at  the  release  location,  but 
other  release  and  atmospheric  conditions  will  result  in  a  mist  pool  that  remains  for  many  minutes  after  the 
tank  empty  time  is  reached. 
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4.0  CONCLUSIONS 


The  mist  pool  theory  predicts  that,  given  appropriate  release  conditions,  the  cloud  formed  from  a  very 
large  release  of  superheated  liquid  over  a  short  release  duration  does  not  move  with  the  wind.  Rather,  the 
vapor  and  liquid  aerosol  collect  at  the  release  location  as  a  low,  two-phase  cloud  and  remain  in  place  on 
flat  terrain.  Vapor  is  entrained  into  the  passing  air  through  the  cloud  top  surface.  The  cloud  mass  is 
progressively  depleted  until  gone.  The  overall  effect  of  formation  of  a  mist  pool  is  that  the  source 
becomes  a  large  area  vapor  flux,  which  lasts  significantly  longer  than  the  two-phase  jet  occurring  while 
the  tank  empties.  Modem  dense  gas  models  do  not  include  methodology  to  simulate  the  processes 
associated  with  mist  pool  formation  and,  as  a  result,  over-predict  concentration  versus  downwind  distance 
for  releases  where  a  mist  pool  forms.  A  source  module  is  needed  that  couples  the  tank  emptying  and  jet 
formation  processes  with  the  mist  pool  entrainment  calculations.  The  module  would  then  provide  an  area 
flux  versus  time  profile  as  input  to  the  transport  and  diffusion  module.  This  paper  provides  a  criterion  to 
predict  whether  a  mist  pool  is  expected  to  form  based  on  release  and  atmospheric  conditions.  The  initial 
methodology  requires  the  user  to  assume  the  air  entrainment  ratio  and  cloud  height,  as  minimal 
supporting  data  exists  to  provide  reliable  estimates.  Rainout  of  droplets  and  reaction/removal  mechanisms 
leading  to  greater  reduction  of  the  cloud  density  have  been  ignored  in  this  analysis. 

It  is  important  to  further  consider  the  mist  pool  concept  because  of  its  dramatic  effect  on  estimation  of  the 
toxic  effects  hazard  area.  Existing  models  predict  that  the  size  of  the  hazard  area  continues  to  increase 
with  release  mass,  supporting  commonly  accepted  observations  about  the  most  effective  atmospheric 
conditions.  Conversely,  the  mist  pool  concept  suggests  that,  for  very  large  releases,  the  toxic  effects 
hazard  area  is  reduced  in  size  due  to  the  extended  duration  release  of  contaminant  at  lower  concentrations. 
The  mass  resulting  in  the  maximum  toxic  effects  area  then  becomes  the  mass  just  less  than  that  which 
results  in  formation  of  a  mist  pool.  The  “most  dangerous”  atmospheric  conditions  of  lower  wind  speed 
and  more  stable  conditions  favor  mist  pool  formation.  As  a  result  of  the  findings  presented  here,  the 
atmospheric  conditions  resulting  in  the  maximum  toxic  effects  area  are  concluded  to  be  shifted  toward 
conditions  of  higher  wind  speed  and  turbulence. 
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6.0  NOMENCLATURE 


Nomenclature 


Symbol  Expanded  Definition 


A  vertical  wind  profile  exponent 

A  cloud  top  area 

B  entrainment  factor 

C  cloud  concentration 

Ci  cloud  concentration  at  the  end  of  stage  1 

Cd  orifice  discharge  coefficient 

Cf  cloud  concentration  after  aerosol  evaporation 

Cp,air  heat  capacity  of  air 

Cpl,h2o  heat  capacity  of  liquid  water 

CpLjic  liquid  heat  capacity 

Cps,H20  heat  capacity  of  ice 

Cpvjic  vapor  heat  capacity 

Cpv,H2o  specific  heat  of  water  vapor 

d  hole  diameter 

Dc  cloud  diameter 

E  volumetric  air  entrainment  ratio 

/  friction  factor 

fc  cloud  height  fraction  when  Richardson  number  equals  12 

FH2o  mass  fraction  of  water  in  ambient  air 

Fj  vapor  mass  fraction  in  jet 

f,  vapor  mass  flux 

Fm,  water  vapor  fraction  remaining  at  final  cloud  temperature 

g  gravitational  acceleration  constant 

Hc  cloud  thickness 

Hd  total  enthalpy  before  depressurization 
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Nomenclature  (cont.) 


Symbol  Expanded  Definition 


Hl 

HLd 

I  f tank 

Htic 

Hv 

Hvd 

L 

Mj 

M2 

Mairid, 

Mc 

Mh20 

Mhs 

Ml 

Mli 

MLd 

Mr 

Mlank 

Mtic 

Mv 


Myc 


MVd 

Mve 

P atm 


liquid  enthalpy  at  boiling  temperature 
liquid  enthalpy  before  depressurization 
total  enthalpy  after  auto-refrigeration 
liquid  enthalpy  at  20  C 
vapor  enthalpy  at  boiling  temperature 
vapor  enthalpy  before  depressurization 
hole  length 

chemical  mass  at  the  end  of  stage  1 

chemical  mass  at  the  end  of  stage  2 

entrained  dry  air  mass 

total  cloud  mass 

entrained  water  mass 

head  space  vapor  mass 

liquid  mass  in  cloud 

chlorine  liquid  mass  at  the  end  of  stage  1 
liquid  mass  in  tank  before  depressurization 
recovered  liquid  mass 
initial  mass 

total  chemical  mass  released 

vapor  mass  in  cloud 

vapor  leaving  the  tank 

vapor  mass  in  tank  before  depressurization 

cloud  aerosol  mass  vaporized 

ambient  air  pressure 
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Nomenclature  (cont.) 


Symbol  Expanded  Definition 


P s.H20 

P  . 

1  s,ice 

P  tank 
P s,  TIC 

P  V 

Q 

Qent 

Qv 

Qvi 

Qv2 

Qv3a 

Qv3b 

Qvf 

R 

Re 

RH 

Ri* 

Ri*f 

t 

t2 

t3a 

t3b 

Tair 

Tb 


saturation  pressure  of  water  at  20  C 

saturation  pressure  of  ice  at  final  cloud  temperature 

tank  pressure 

chemical  saturation  pressure 

partial  pressure  of  chlorine  vapor  in  cloud 

jet  mass  flow  rate 

average  entrainment  rate 

mass  entrainment  rate 

average  mass  entrainment  rate  for  stage  1 

average  mass  entrainment  rate  for  stage  2 

average  mass  entrainment  rate  for  stage  3,  part  1 

average  mass  entrainment  rate  for  stage  3,  part  2 

mass  entrainment  rate  after  aerosol  evaporation 

ideal  gas  constant 

jet  Reynold's  number 

relative  humidity 

Richardson  number 

Richardson  number  after  aerosol  evaporation 

time  from  hole  formation 

duration  of  stage  2 

duration  of  stage  3,  part  1 

duration  of  stage  3,  part  2 

ambient  air  temperature 

chemical  boiling  temperature 
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Nomenclature  (cont.) 


Symbol  Expanded  Definition 


Tc 

tent 

tf 

Tfmo 

Ttank 

U* 

V 

Vair 

Vc 

Vh20 

Vice 

VL 

V tank 

Vtic 

Vv 

w* 

w  ■ 

VY  air 

Wh20 

W TIC 

a 

AHfmo 

Ally.mo 

AH v,  tank 

Ally  TIC 


final  cloud  temperature 
entrainment  duration 
tank  empty  time 
freezing  temperature  of  water 
tank  temperature 
ambient  air  friction  velocity 
jet  velocity 

volume  of  air  entrained 
cloud  volume 
water  vapor  volume 
ice  volume 

chemical  liquid  volume 
tank  volume 

volume  of  chemical  mass  released 
chemical  vapor  volume 
vertical  entrainment  velocity 
molecular  weight  of  air 
molecular  weight  of  water 
chemical  molecular  weight 
mass  evaporation  ratio 
heat  of  fusion  of  water 
heat  of  vaporization  of  water 
heat  of  vaporization  at  20  °C 

chemical  heat  of  vaporization  at  boiling  temperature 


24 


Nomenclature  (cont.) 


Symbol 

Expanded  Definition 

AP 

pressure  drop  along  hole  length 

AT 

temperature  decrease  in  tank 

AU air, dry 

energy  change  from  cooling  air 

AUh20 

energy  change  from  cooling  water  vapor 

AU  ice 

energy  change  from  forming  and  cooling  ice 

AUl 

energy  change  from  cooling  chemical  liquid 

AUV 

energy  change  from  cooling  chemical  vapor 

Mtic 

liquid  viscosity  of  chlorine 

Pair 

ambient  air  density 

Pair, dry 

partial  density  of  dry  air 

Pc 

cloud  density  with  entrained  air 

Pci 

cloud  density  at  the  end  of  stage  1 

pcf 

cloud  density  after  aerosol  evaporation 

pcL 

liquid  density  at  final  cloud  temperature 

PH20 

partial  density  of  water  vapor 

Pice 

density  of  ice 

PJ 

density  of  jet  at  boiling  temperature 

PL 

liquid  density  at  boiling  temperature 

PLd 

liquid  density  depressurization 

PL,tank 

liquid  density  in  the  tank 

Ptank 

fluid  density  in  the  tank 

Pv 

vapor  density  at  boiling  temperature 

PVd 

vapor  density  before  depressurization 
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